Critical Behavior of Three-Dimensional Disordered Potts Models with Many States 
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We study the 3D Disordered Potts Model with p = 5 and p — 6. Our numerical simulations (that 
severely slow down for increasing p) detect a very clear spin glass phase transition. We evaluate 
the critical exponents and the critical value of the temperature, and we use known results at lower 
p values to discuss how they evolve for increasing p. We do not find any sign of the presence of a 
transition to a ferromagnetic regime. 



PACS numbers: 75.50.Lk, 75.40.Mg, 64.60.F- 
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I. INTRODUCTION 

The three dimensional (3D) disordered Potts model 
(DPM) is an important system, that could help in clari- 
fying a number of open and crucial questions. The first 
issue that comes to the mind is the possibility of under- 
standing the glass transition, since this is a very chal- 
lenging problem. On more general grounds, it is very 
interesting to try and qualify the behavior of the sys- 
tem when the number of states p becomes large: here 
we should see the paradigm of a "hard" , first order like 
transition but, as we will discuss in the following, only 
sometimes this turns out to be clear (see for example the 
set of large scale, very accurate numerical simulations of 
Ref. |l|, dealing with a model slightly different from the 
one defined here). 

In such a difficult situation extensive numerical sim- 
ulations are more than welcome, and the Janus super- 
computer— i 3 -, optimized for studying spin glasses, reaches 
its peak performances when analyzing lattice regular sys- 
tems based on variables that can take a finite, small num- 
ber of values: disordered Potts models fit very well these 
requirements. Using the computational power of Janus 
we have been able to consistently thermalize the DPM 
with p = 5 and 6 on 3D (simple cubic) lattice systems 
with periodic boundary conditions and size up to L = 12. 
Bringing these systems to thermal equilibrium becomes 
increasingly harder with increasing number of states: it 
has been impossible for us, even by using a large amount 
of time of Janus (that for these problems performs, as we 
discuss better in the following, as thousands of PC pro- 
cessors), to get a significant, unbiased number of samples 



thermalized, and reliable measurements of physical quan- 
tities, for p > 5 on a L = 16 lattice. 

Our results lead us to the claim that the critical be- 
havior of the DPM with a large number of states p is very 
subtle, and if p is larger than, say, 5, numerical simula- 
tions could easily give misleading hints. The numerical 
results that we will discuss in the following lead us to 
believe that the spin glass transition gets stronger with 
increasing number of states p: a theoretical analysis of 
these results suggests that the transition could eventually 
become of first order for p large enough. We do not ob- 
serve, for both p — 5 and p = 6, any sign of the presence 
of a spontaneous magnetization. 



II. MODEL AND OBSERVABLES 

We have performed numerical simulations of the DPM 
on a simple cubic lattice of linear size L with periodic 
boundary conditions. The Hamiltonian of the DPM is 
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where the sum is taken over all pairs of first neighboring 
sites. In the p-states model spins Si can take p different 
values {0, 1, . . . , p— 1}. In this work we analyze the p — 5 
and 6 cases. The couplings Jij are independent random 
variables taken from a bimodal probability distribution 
(Jij = ±1 with probability ^). For a different definition 
of a disordered Potts model see Ref. |3- 

It is convenient to rewrite the variables of the Potts 
model using the simplex representation, where the p Potts 
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states are described as vectors pointing to the corners 
of a (p — 1) dimensional hyper-tetrahedron. The Potts 
scalar spins Sj are thus written as (p — l)-dimensional 
unit vectors Si satisfying the relations 
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where a and b € [l,p]- We use this vector representation 
to define the observables required to investigate the criti- 
cal behavior of the system. In the simplex representation 
we have that: 
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The couplings in the simplex representation have the 
form 
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P-1 
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The spin glass behavior is studied via a properly de- 
fined tensorial overlap between two replicas (independent 
copies of the system characterized by the same quenched 
disorder variables Jij). Its Fourier transform (with wave 
vector k) is given by£ 

«""(*) 
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where S^ 1 ^ is the /i component of the spin at site i of 
the first replica in the simplex representation, S^ 1 ' the 
v component of the spin at site i in the second replica, 
and V — L 3, is the volume of the system. 

This spin glass order parameter is then used to define 
the spin glass susceptibility in Fourier space. 
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where ((• • •)) indicates a thermal average and (• • •) de- 
notes the average over different realizations of the disor- 
der (samples in the following). With the above definition, 
X ? (0) is the usual spin glass susceptibility. 

We are interested in studying the value of the dimen- 
sionless correlation length £/£, since at the transition 
temperature it does not depend on L, and is therefore ex- 
tremely helpful to estimate the critical temperature value 
T c : in fact one can usually simulate different lattice sizes, 
and look for the crossing point in the plot of the different 
£/L values. One can derive^ the value of the correlation 
length £ from the Fourier transforms of the susceptibility 
with 
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where fcm is the minimum wave vector allowed in the 
lattice. With the periodic boundary conditions used in 



this work we have fem = (2ir/ L,0,0) or any of the two 
vectors obtained permuting the indexes. 

We also study the ferromagnetic properties of the 
model by monitoring the usual magnetization 



m 



and correspondingly the magnetic susceptibility 

Xm = v{K). 



(8) 
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These two observables are crucial to check the possible 
existence of a ferromagnetic phase, as predicted by the 
mean field approximation of this model^. 



III. NUMERICAL METHODS 

We have analyzed the DPM with 5 and with 6 states, 
on a number of lattice sizes (L = 4, 6, 8, and 12). All the 
numerical simulations have been run using a standard 
Metropolis algorithm combined with the Parallel Tem- 
pering (PT) optimized algorithm, in order to improve 
performances and allow to reach thermalization despite 
the very large relaxation times typical of spin glass mod- 
els. 

We define a Monte Carlo sweep (MCS) as a set of 
V trial updates of lattice spins. Each simulation con- 
sists on a thermalization phase, during which the sys- 
tem is brought to equilibrium, and a phase of equilib- 
rium dynamics in which relevant physical observables are 
measured. As we require high quality random numbers, 
we use a 32-bit Parisi-Rapuano shift register— pseudo- 
random number generator^ 

In order to improve the simulation performance and 
to speed up thermalization we apply a step of the PT al- 
gorithm^ every few MCS's of the Metropolis algorithm. 
The PT algorithm is based on the parallel simulation of 
various copies of the system, that are governed by differ- 
ent values of temperature, and on the exchange of their 
temperatures according to the algorithm's rules. In prac- 
tice we let the different configurations evolve indepen- 
dently for a few MCS, and then we attempt a tempera- 
ture swap between all pairs of neighboring temperatures: 
the aim is to let each configuration wander in the allowed 
temperature range (that goes from low T values, smaller 
than T c , to high T values, larger than T c ), and to use the 
decorrelation due to the high T part of the landscape to 
achieve a substantial speed up. 

In order to check the time scales of the dynamical pro- 
cess, so as to assess the thermalization and the statistical 
significance of our statistical samples, we have computed 
a number of dynamical observables that characterize the 
PT dynamics. 

One of them is the temperature-temperature time cor- 
relation function, introduced in Ref. 10, that we briefly 
recall. Let /3^(t) be the inverse temperature of the 
system i at time t (i = 0, ...,Nt — 1), where Nt is 
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the total number of systems evolving in parallel in the 
PTi^i We consider an arbitrary function of the sys- 
tem temperature, /(/?), changing sign at f3 c . We shall 
name / t = f(/3^(t)). In equilibrium, system i can 
be found at any of the Nt with uniform probability, 
hence (/ t W ) = Ef-To" 1 f{Pk)/N T , for all % and all t. We 
must choose a function / as simple as possible, such that 
X)fc=o 1 fifth) = Next, we can define the correlation 
functions 



cf{t) = 
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where N is the total simulation time. To gain statistics 
we consider the sum over all the systems 



JVt-1 
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Notice that this correlation function measures correla- 
tions for a given copy of the system, that is characterized, 
during the dynamics, by different temperature values. 

We have characterized the correlation function Pf(t) 
through its integrated autocorrelation time&ii: 



Tint 



dt p f (t) , 
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where A; nt = w r; n t and we have used u — 10 (we have 
always used a total simulation time larger than 15 or 20 
times Ti n t). 

We have studied the systems defined on the smaller 
lattices (L = 4 and 6) on standard PCs, while for the 
larger lattice sizes we have used the Janus computer— 
an FPGA-based machine specifically designed to handle 
simulations of spin glass models. The performance im- 
provement offered by Janus allowed us to thermalize lat- 
tices of size up to L = 12. While the thermalization 
of lattices with L = 8 was relatively fast, the bigger lat- 
tice sizes proved to be rather difficult to equilibrate, even 
within Janus, things getting worse as the number of Potts 
states increases. 

Tables Q] and [TT] summarize the details about the nu- 
merical simulations respectively for the p = 5 and the 
p = 6 case. We were able to thermalize a large number 
of samples for L up to 12. The thermalization of L = 16 
is possible, but it requires a dramatically large invest- 
ment in computer resources, since the time required by 
each sample is very large. Because of that, and given the 
resources we could count upon, we have only been able 
to analyze a few samples: the results for the few sam- 
ples that we have studied in this case are consistent with 
the ones obtained from the smaller sizes. In addition, 



L 


■^Vj samples 


MCS mi n 


[/3min 7 /^max] 


Np 


-^Metropolis 


Nm 


4 


2400 


10 7 


[1.6, 9.5] 


18 


5 


10 3 


6 


2400 


2 x 10 7 


[1.6, 9.5] 


22 


5 


10 3 


8 


2448 


4 x 10 8 


[1.7, 6.5] 


24 


10 


2 x 10 5 


12 


2451 


6 x 10 9 


[1.8, 5.5] 


20 


10 


2 x 10 5 



TABLE I: Details of the simulations for p — 5. N 3a .mples is 
the number of samples (i.e. of the disorder realizations that 
we have analyzed), MCSmin is the minimum number of MCSs 
that we have performed, [/3 m in, /3max] is the range of inverse 
temperatures simulated in the PT, Np is the number of tem- 
peratures inside this interval, -^Metropolis is the frequency of 
the Metropolis sweeps per PT step, and iVm is the total num- 
ber of measurements performed within each sample. 



L 


N samples 


MCS min 


[/3miii7 /^max] 


N p 


-^Metropolis 


N m 


4 


2400 


10 7 


[2.1, 9.8] 


10 


5 


10 3 


6 


2400 


2 x 10 7 


[2.0, 9.65] 


16 


5 


10 3 


8 


1280 


10 9 


[1.7, 7.5] 


30 


10 


2 x 10 5 


12 


1196 


6 x 10 10 


[1.6, 6.5] 


22 


10 


2 x 10 5 



TABLE II: As in table □ but for p = 6. 



for some samples with L — 8 and L — 12, which were 
especially difficult to thermalize, we had to use larger 
numbers of MCS's: see section HV AI 

The number of Metropolis sweeps per PT step is 10 on 
Janus and 5 on the PC, and there is an important reason 
for that: in a standard computer the time needed for a 
step of the PT algorithm is small compared with a com- 
plete Metropolis MCS. This is not true on Janus, where 
it takes longer to perform a PT step than an Metropolis 
MCS: because of that, after a careful test of the over- 
all simulation performance, we decided to lower the PT 
to Metropolis MCS ratio in order to increase Janus effi- 
ciency. 

In the p = 5 case a numerical simulation of a sin- 
gle sample (thermalization plus measurements) on Janus 
takes 39 minutes for L = 8 and 10 hours on L = 12. The 
same simulations would require 7.4 days of an Intel^ 
Core2Duo ( ™> 2.4 GHz processor for L = 8 and 315 
days for L — 12. These values grow when p = 6: here 
the equilibration takes 120 minutes for an L — 8 sample 
and 110 hours for L = 12 (on the PC they would take 24 
days for L — 8 and 10 years for L = 12). 

The results shown in this paper for the p = 5 model 
would have required approximately 2150 equivalent years 
ofanIntel( fl ) Core2Duo(™) 2.4 GHz processor: the ones 
for p = 6 would have required 12000 years. 
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IV. RESULTS 
A. Thermalization Tests 

Thermalization tests are a crucial component of spin 
glass simulations. Before starting to collect relevant re- 
sults from the data we have to be sure that they are 
actually taken from a properly thermalized system, and 
are not biased from spurious effects. 

A standard analysis scheme consists in evaluating the 
average value of an observable on geometrically increas- 
ing time intervals. The whole set of measurements is 
divided in subsets, each of which covers only part of the 
system's history (the last bin covers the last half of the 
measurements, the previous bin takes the preceding quar- 
ter, the previous bin the previous eighth and so on), and 
observables are averaged within each bin. The conver- 
gence to equilibrium is checked comparing the results 
over different bins: stability in the last three bins within 
error bars (that need to be estimated in an accurate way) 
is a good indicator of thermalization. 

We show in figures [1] and [5] the logarithmic binning of 
£, as defined in equation , in the p = 5 and p — 6 cases. 
The compatible (and stable) values for the three last 
points satisfy the thermalization test explained above. 
The data in the plots are for the lowest temperature used 
on each lattice size: this is expected to be the slowest 
mode of the system, and its thermalization guarantees 
that also data at higher temperature values are thermal- 
ized. The plateau in the last part of each plot is a clear 
signal of proper thermalization: only data from the last 
bin are eventually used to compute thermal averages. 
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FIG. 1: Log-binning thermalization test for p — 5. For all 
data points the point size is bigger than the corresponding 
error bar. 

We have also investigated how thermalization is 
reached in the individual samples (as opposed to the in- 
formation on averages obtained from figures [I] and [5]) : 
to do that we have studied the correlation function for 
the temperature random walk defined in (|12[) and its as- 
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FIG. 2: As in figure [2 but p = 6. 
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FIG. 3: The autocorrelation function (|12[) for one generic 
sample (p = 6, L = 8). 

sociatcd integrated autocorrelation time, T- lnt , defined in 
(fT5)) . As an example we plot in figure [3] the autocorrela- 
tion function (fT2)) for a given sample as a function of the 
Monte Carlo time (here L = 8 and p = 6): one can see a 
fast, exponential decay in the left part of the figure, and 
(large) fluctuations around zero at later times. 

Sample to sample fluctuations of 7i n t are very large: in 
figure [4] we plot T lnt for all our samples with p — 5, L = 8. 
In order to be on the safe side we have increased the 
number of MCS, by continuing the numerical simulation 
for a further extent, in all samples where our estimate of 
Tj nt was bigger than the length of the simulation divided 
by a constant c (c = 20 for L = 8 and c = 15 for L = 12, 
where achieving thermalization is much more difficult). 23 

B. Critical temperature and critical exponents 

Our analysis of the critical exponents of the system 
has been based on the quotient method^'i^: by using the 
averaged value of a given observable O measured in lat- 
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FIG. 4: Integrated autocorrelation time, Tint, for all p = 5, 
L = 8 samples. r m t is in units of blocks of ten measurements, 
i.e. of 2010 3 MCS. Samples above the green line have been 
"extended" (see the text for a discussion of this issue). 



FIG. 5: Overlap correlation length in lattice size units as a 
function of the inverse temperature /3 for L = 4, 6, 8 and 12. 
Here p = 5. 



tices of different sizes, we can estimate its leading critical 
exponent xq, 



(0(/3)) « |/3 - & 



(14) 



By considering two systems on lattices of linear sizes L 
and sL respectively one has that&i^ 



m,sL)) 



= s xo,w 



+ 0(L-").. 



(15) 



where v is the critical exponent of the correlation length 
and uj is the exponent of the leading-order scaling- 
corrections^. 

We use the operators c^£, from ([7]), and x 9 , from © 
in equation (|15[) to obtain respectively the critical ex- 
ponents l + l/i/ and 2 — rj q . The exponent 2 — 7y m is 
obtained applying eq. (|15j) to the magnetic susceptibility 
Xmi from ©. 

To use the quotient method we start estimating the 
finite-size transition temperature: we do this by looking 
at the crossing points of the correlation length in lattice 
units (£/£) for various lattice sizes. We have used a 
cubic spline interpolating procedure to compute both the 
crossings of £/L and its /3-derivative (we have followed 
the approach described in detail in Ref. [l3h . 

We show in figures [5] and [6] the behavior of £/£ as a 
function of /3. The different curves are for different lat- 
tice sizes. The crossing points are rather clear in both 
cases, giving a strong hint of the occurrence of a second 
order phase transition. At least for p — 5 scaling correc- 
tions play a visible role, and the crossing points undergo 
a small but clear drift towards lower temperatures for in- 
creasing lattice sizes. We summarize in tables Hill and ITVl 
the (3 values of the crossing points for two different pairs 
of lattice sizes, together with the estimated values of the 
critical exponents v and r\ q that we obtain using relation 




FIG. 6: As in figure El but p = 6. 



Since we can only get reliable results on small and 
medium size lattice we cannot control in full scaling cor- 
rections, and a systematic extrapolation to the infinite 
volume limit is impossible. It is clear however that the ef- 
fective critical exponents summarized in tables lnil and irvl 
do not suggest that asymptotically for large volume the 
system will not be critical (in this case, for example, 
•q q should be asymptotically equal to 2): our numerical 
data clearly support the existence of a finite temperature 
phase transition. 





Across (Ll , L2) 


u(L-i,L 2 ) 


ri q (Li,L 2 ) 


■rj m (Li, L 2 ) 


(4,8) 


4.83(5) 


0.82(3) 


0.13(2) 


1.72(2) 


(6,12) 


5.01(4) 


0.81(2) 


0.16(2) 


1.94(2) 



TABLE III: Numerical values of our estimates for the crossing 
point of the curves £/L. We give Across, the thermal critical 
exponent v, the anomalous dimension of the overlap rj q , and 
the anomalous dimension of the magnetization r\ m . 
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(Li,L,2, 
(4,8) 
(6, 12) 



/3cross(Jjl, L2) 

6.30(9) 
6.26(7) 



u(Li,L2) 
0.80(2) 
0.80(4) 



T}q (1/1,1/2) 

0.10(2) 
0.16(2) 



1.453(19) 
1.971(19) 



TABLE IV: As in table EH but p = 6. 

We take as our best estimates for the critical exponents 
the one obtained from the lattices with sizes L = 6 and 
L = 12. For p = 5 

/3 C = 5.01(4) , v = 0.81(2) , Vq = 0.16(2) , (16) 

while for p = 6. 

& = 6.26(7) , 1/ = 0.80(4) , ^ = 0.16(2) . (17) 

It is interesting to compare these values with those of 
other Potts models with a different number of states. In 
particular we are interested in the value of the critical 
exponents as a function of the number of states, since we 
want to characterize the critical behavior of the various 
models and attempt a prediction of the model's behav- 
ior when the number of states is large. In our particular 
model and with the (low) values of the temperature that 
are interesting for us (since we need to get below the 
critical point) even with the large computational power 
available to us thanks to Janus the simulation for p = 8, 
say, on a L = 12 lattice, would require an unavailable 
amount of CPU time. What is found in the very inter- 
esting work of Refs. [l] and[H is different, since there one is 
able to thcrmalize a p = 10 model on a large lattice, and 
no transition is observed. The model analyzed in these 
two references^ is indeed slightly (or maybe, it will turn 
out, not so slightly) different from the present one, since 
there J is negative. It is not clear to us if this difference 
could explain a quite dramatic discrepancy of the ob- 
served behavior, or if, for example, a different (very low) 
temperature regime should be analyzed to observe rele- 
vant phenomena: this is surely an interesting question 
to clarify, and the fact that the coupling have a negative 
expectation value, reducing in this way frustration, could 
turn out to make a difference. 



C. Absence of ferromagnetic ordering in the 
critical region 

Our DPM is in principle allowed to undergo a ferro- 
magnetic phase transition (since no symmetry protects 
it) , and at low temperatures could present a spontaneous 
magnetization, as discussed in Ref. [l3| . Because of that 
we have carefully studied the magnetic behavior of the 
model at low temperatures. We have analyzed both the 
magnetization and the magnetic susceptibility below the 
spin glass critical point. 

In the paramagnetic phase the magnetization is ran- 
dom in sign, and its absolute value is expected to be pro- 
portional to 1/y/V. In Figs. [7] and [5] we check whether 
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( T - 1 

123456789 10 

P 

FIG. 7: Magnetic susceptibility as a function of /3 for L = 4, 
6, 8 and 12. Here p = 5. 
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FIG. 8: As in figure [3 but p = 6. 



(|m|) around the spin glass critical region tends to an 
asymptotic value for larger lattice size, or not. From 
the figures we see (|m|) goes to zero in the critical 
region. Also, we studied the magnetic susceptibility 
Xm = y(\ m \' 2 ) which is independent of size. Again in 
Figs. [7] and [8] we check that, and we see a non-divergent 
behavior. This behavior is extremely different from a 
ferromagnetic phase in which \m diverges as the volume. 

Besides, as reported in Sec. IIVBI the exponent rj m 
is close to 2, so we could say that a ferromagnetic- 
paramagnetic phase transition does not happen in the 
range of temperatures that we have studied. 



V. EVOLUTION OF CRITICAL EXPONENTS 
WITH p 

In table [V] we summarize the values of the inverse crit- 
ical temperature and of the thermal and overlap critical 
exponents for DPM from p = 2 (the Ising, Edwards- 
Anderson spin glass) up to p — 6. We also plot these 
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data items in figure [9l 

From table [V] and figure |H] some results emerge very 
clearly. First, the inverse critical temperature roughly 
follows a linear behavior in p, with a slope very close to 
one. We have added in table|V]the ratio (R) between the 
numerical determinations (in 3d) of f3 c {p) and their values 
in the Mean Field (MF) approximation. One can see that 
the large deviations from the MF prediction occur for 
large values of p (notice that R > 1 since MF suppresses 
fluctuations) .— 

Second, v decreases monotonically and r\ q grows mono- 
tonically with the number of states p. To discuss this be- 
havior it is useful to keep in mind that when using finite 
size scaling to study a disordered first order phase tran- 
sition one expects to findii v — 2/D and 2 — r/ q = D/2, 
i.e., in our D = 3 case, v = 2/3 and r\ q — 1/2. These 
are "effective" exponents, that are a bound to the ones 
allowed for second order phase transitions. 

Both sets of values for v and r\ q are indeed completely 
compatible with tending, as p increases, to those limit 
values that characterize a first order phase transition. If 
this turns out, as our numerical data make very plau- 
sible to be true, two different scenarios open. The first 
possibility is that the p-states DPM undergoes a disor- 
dered first order phase transition for large enough values 
of p (just as in the ordered Potts model, that for p > 3 
undergoes a first order phase transition), while the sec- 
ond possibility is that the DPM will show a standard 
second order phase transition for all finite values of p. 
This is the typical issue that is very difficult to settle 
with numerical work: an analytical solution of the model 
with infinite number of states would be very useful as a 
starting point in order to discriminate between these two 
possible scenarios. 



p 


/3c 


V 


Vq 


R 


2 (Ref.fl5]) 


1.786(6) 


2.39(5) a 


-0.366(16) 6 


2.187(8) 


2 (Ref.[16]) 


1.804(16) 


2.45(15) 


-0.375(10) 


2.209(20) 


3 (Ref.[5j) 


2.653(35) 


0.91(2) 


0.02(2) 


2.17(3) 


4 (Ref.[13j) 


4.000(48) 


0.96(8) 


0.12(6) 


2.45(3) 


5 (this paper) 


5.010(40) 


0.81(2) 


0.16(2) 


2.51(2) 


6 (this paper) 


6.262(71) 


0.80(4) 


0.16(2) 


2.69(3) 



"This value of v is from It is different and more reliable than 

the one obtained from the spin glass susceptibility, that, because 
of large scaling corrections, would severely depend on the kind of 
analysis. 

'This value of r\ q is from the study of the spin glass susceptibility. 

TABLE V: Critical parameters as a function of p. All data 
are for binary couplings, with zero expectation value. By R 
we denote the ratio between the critical /3 in three dimensions 
and that computed in Mean Field. 
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FIG. 9: In the bottom plot: f3 c versus p, and the straight line 
f(p) — p. Middle plot: v as a function of p. We also show 
(dashed line) the value which marks the onset of a disordered 
first order phase transition (z/fi rs t = 2/3). Upper plot: r\ q as a 
function of p. 



VI. CONCLUSIONS 

In this note we have characterized the critical behav- 
ior of the 3D DPM with p = 5 and p = 6, i.e. with a 
reasonably large number of states. Our numerical sim- 
ulations have allowed us to reach some clear evidences, 
and to stress some difficult issues that will require further 
analysis. 

We first stress that in both cases the spin glass tran- 
sition is very clear, and we have been able to obtain a 
reliable estimate of the critical temperature and of the 
critical exponents v and rj q . We have discussed what 
happens when p increases; we have found that /3 C in- 
creases like p. A similar result was conjectured in Ref. 
[H (for all values of p) analyzing high temperature series 
and found in Mean Field for p < 4 (although, of course, 
the slope is wrong). In addition, the behavior of v and r\ q 
is compatible with going to the large p limit value that 
characterizes a first order phase transition. 

In the low temperature regime we do not see any sign of 
a transition to a ferromagnetic regime, that would be in 
principle allowed by the structure of our model. We can- 
not exclude that at very low T values something would 
happen, but in all the range we can explore the system 
stays in the spin glass phase. 

A last piece of important evidence is that low temper- 
ature simulations of this model look difficult, and that 
they slow down severely for increasing p. In our particu- 
lar model, where the expectation of the coupling is zero, 
it would be impossible to study reliably a p = 8 model 
with the computational resources available today. 

This last observations opens indeed a last point that 
it will be interesting to analyze in the future. When cou- 
plings have a negative expectation value the simulation 
of a p = 10 model— £ is possibly easier than it would be in 



8 



our case, and the results are very different: in that case 
one does not see any sign of a phase transition. Analyz- 
ing how the DPM depends on the expectation value of 
the couplings is indeed at this point a crucial issue, since 
it could turn out that the reduction in frustration due 
to a negative net value of the couplings could completely 
change the critical behavior of the model. 
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that T a /J = 1 for p < 4 and (T c /J) 2 = 1 + (p - 4) 2 /42 + 
0((p — 4) 4 ) for p > 4. In addition for very large p, T c /J — 
i (p/ logp) 1 / 2 . Taking into account the extra p factor in 
the Hamiltonian used in the Mean Field and the fact that 
J = \pM (Jfj = J 2 /N, being N the number of spins in the 
MF computation) since we are working in finite dimension 
(d), we obtain the finite dimension version of the critical 
P using the Mean Field approximation: f3 c — p/v2d for 
p < 4 and p e = -fa (1- (p-4) 2 /84 + 0((p-4) 4 )) for 

p > 4 (notice the minus signum of the (p — 4) 2 correction); 

in addition, for large p, one obtains /3 C ~ (plogp) 1//2 . 

Note that in our case %/2d ~ 2.45. 



